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Abstract 

We present a Lorentz-invariant formulation of baryon chiral perturbation theory including spin-3/2 fields. 
Particular attention is paid to the projection on the spin-3/2 components of the delta fields. We also discuss 
the nucleon mass and the pion-nucleon sigma term. 
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1. The delta resonance plays a special role in low-energy nuclear and particle physics, due to its near mass 
degeneracy with the nucleon and its strong couplings to pions, nucleons and photons. It was therefore argued 
early that spin-3/2 (decuplet) states should be included in baryon chiral perturbation theory p^, which is the 
low-energy effective field theory of the Standard Model. The work of 1 and subsequent authors made use of the 
heavy baryon approach, which treats the baryons as static sources like in heavy quark effective field theory and 
allows for a systematic power counting in the presence of matter fields, as pioneered in [2] and systematically 
explored in [3] . Also, special care has to be taken about the decoupling of resonances in the chiral limit 4 . This 
approach was systematized by counting the nucleon-delta mass splitting as an additional small parameter in 
Ref. |3] , the corresponding power counting was called the "small scale expansion" . The heavy baryon approach 
has been successfully applied to a variety of processes, for reviews see d, and a status report on chiral 
effective field theories with deltas is given in More recently, a Lorentz-invariant formulation of baryon 
chiral perturbation theory has become available 9 , the so-called "infrared regularization" (IR). A similar 
approach had been presented earlier in |10|. Such a formulation allows to include strictures from analyticity 
and is thus particularly suited for extensions of chiral perturbation theory based on dispersion relations, see 
e.g. El E21 C31 EI EH EH El- Furthermore, the use of the Dirac propagator as opposed to the static fermion 
propagator in the heavy baryon scheme allows for a resummation of important recoil effects, see e.g. the work 
on electromagnetic form factors JS] or the nucleon spin structure ^3 EH]- The chiral expansion of the baryon 
masses also seems to converge faster in this scheme [21] whereas the description of pion-nucleon scattering is 
not yet in a satisfactory status [221 123| . The renormalization of relativistic baryon chiral perturbation theory 
has been discussed in detail in 24 . In this Letter, we give a consistent extension of the infrared regularization 
method in the presence of spin-3/2 # 4 . What is new here is that we explicitly project onto the spin-3/2 parts 
and that the method easily allows to include any external source in a chiral and gauge invariant fashion (in 
contrast to the recent proposal in |2H]). As simple applications, we calculate the delta loop contribution to the 
nucleon self-energy and the pion-nucleon sigma term. 

2. In this section, we briefly review the formalism necessary to describe spin-3/2 fields. More details can e.g. 
be found in |2fill27ll2"%| . We first write down the propagator for a free spin-3/2 field (called the delta propagator 
from here on) in d space-time dimensions, 

A i \ rf+mA J 1 (d-2)p ll p v - p„7 M 1 

G ^ {P) = -p^ni -J—l 1 ^ - {d-l) m \ + (d-l)m A ) ' (1) 

with toa the delta mass. This may be rewritten by projection onto the spin 3/2 and 1/2 components 

w - ^ - vi^: + « 2 ),,) + j^h ^ ~> W\„ ■ w 

with 

d3/2 1 1 t u , a d-Ap^ 

P„L = 9^ - — J lulv - (d _ 1)p2 [tflnPv + PulvP) - J—£ — , 

These spin projection operators fulfill the orthogonality relations 

{P^^PitY" = 8 IJ 5 jk {P( l ) P ix , (I, J) = 1,2, (i,j,k,l) = 1,2. (4) 

Note the infrared singular pieces ~ 1/p 2 appearing in the spin-projected parts of the propagator, which will 
play an important role later on. Also, we remark that the spin-1/2 pieces do not propagate and thus one should 

# 4 The inclusion of the delta when one separates loop integrals into soft and hard parts was already considered in 1231 . 
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be able to absorb their contribution in purely polynomial terms (which amounts to a redefinition of certain 
low-energy constants in the effective field theory). If one considers processes where the center-of-mass energy 
reaches the delta mass, \/s — m&, the free delta propagator has to be resummed as described in detail e.g. in 
^3 ESI- Let us now briefly discuss the leading chiral pion-delta-nucleon Lagrangian. It is given by 

Can = c A % 9£ (Z)< ^ + h.c. , (5) 

where ^ describes the Rarita-Schwinger (spin-3/2) field, ip the nucleon doublet, = Tr(Tjtt^)/2 the axial 
current and O^^ = + (Z — 1/2)7^7,,. The dependence on the so-called off-shell parameter Z is, however, 
spurious, since it does not appear in the spin-3/2 contributions and thus physical observables do not depend on 
it. The axial- vector coupling constant ca, which is related to the pion-delta-nucleon coupling g n AN, is frequently 
called h,A in the literature. More precisely, what appears in Eq.JSJ) is the the coupling ca in the chiral limit. 
For simplicity, we use in the following the SU(4) coupling constant relation 

3 

c A = h A = 9A ■ (6) 

The resulting value is somewhat larger than found in typical fits say to pion-nucleon scattering data or from 
fitting the A — > Nir width, see e.g. jSHHSIIj- We remark that like in Eq.JSJ), we are using standard chiral-invariant 
couplings of the it AN system, because differences to other couplings appearing in the literature can simply be 
absorbed in the polynomial contributions to be discussed in the following. 



3. Next, we make some short remarks about the effective field theory of massive spin- 1/2 fields chirally coupled 
to Goldstone bosons and external sources, called baryon chiral perturbation theory. It is complicated by the 
fact that the nucleon mass does not vanish in the chiral limit and thus introduces a new mass scale apart from 
the ones set by the quark masses. Therefore, any power of the quark masses can be generated by chiral loops 
in the nucleon (baryon) case, spoiling the one-to-one correspondence between the loop expansion and the one 
in the small parameter q. One method to overcome this is the heavy mass expansion (called heavy baryon 
chiral perturbation theory, for short HBCHPT) where the nucleon mass is transformed from the propagator 
into a string of vertices with increasing powers of 1/m. Then, a consistent power counting emerges^ 5 . However, 
this method has the disadvantage that certain types of diagrams are at odds with strictures from analyticity. 
The best example is the so-called triangle graph, which enters e.g. the scalar form factor or the isovector 
electromagnetic form factors of the nucleon. In a fully relativistic treatment, such constraints from analyticity 
are automatically fulfilled. It was argued in that relativistic one-loop integrals can be separated into "soft" 
and "hard" parts. While for the former the power counting as in HBCHPT applies, the contributions from 
the latter can be absorbed in certain LECs. In this way, one can combine the advantages of both methods. A 
more formal and rigorous implementation of such a program was given in [55]. The underlying method is called 
"infrared regularization" . Any dimensionally regularized one-loop integral H is split into an infrared singular 
(called /) and a regular part (called R) by a particular choice of Feynman parameterization, 

H = I + R. (7) 

Consider first the regular part. If one chirally expands the contributions to R, one generates polynomials in 
momenta and quark masses. Consequently to any order, R can be absorbed in the LECs of the effective 
Lagrangian. On the other hand, the infrared (IR) singular part has the same analytical properties as the full 
integral H in the low-energy region and its chiral expansion leads to the non-trivial momentum and quark-mass 
dependences of CHPT, like e.g. the chiral logs or fractional powers of the quark masses. For a typical one- 
loop integral (like e.g. the nucleon self-energy £) this splitting can be achieved in the following way (omitting 
prefactors) 

d d k 1 f 1 , f d d k 1 

dz 



(2n) d AB J J (27r) d [(1 - z)A + zBf 

d d k 



OO />OC 



h J J (2n) d [(l-z)A + zBf 



= I + R, (8) 



# 5 The extension of this method for including deltas treating the NA mass splitting as an additional small parameter has been 
given in 0. 
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with A = M 2 ~ k 2 — ie, B = m 2 — (p — k) 2 — ie, e — ► + , M T the pion mass and d the number of space-time 
dimensions. Any general one-loop diagram with arbitrary many insertions from external sources can be brought 
into this form by combining the propagators to a single pion and a single nucleon propagator. It was also shown 
that this procedure leads to a unique, i.e. process-independent result, in accordance with the chiral Ward 
identities of QCD j)- . Consequently, the transition from any one-loop graph H to its IR singular piece I defines 
a symmetry-preserving regularization. For more details, the reader is referred to [Hj • 

4. Next, we develop a systematic infrared regularization for effective field theories with spin-3/2 fields. This 
should be considered the main result of this paper. Due to the presence of the IR singular terms ~ 1/p 2 when 
one projects the delta propagator onto its spin-3/2 and spin-1/2 components, cf. Eq.©, one has to deal with 
a type of integrals that do not appear in the pure pion-nucleon approach. To be specific, consider the integral 
appearing in the self-energy (details of the self-energy calculation will be given below) 



J°o(p 2 



i r d^k -i 



i J (2tt) 4 {k-p) 2 {M 2 -k 2 ) 
dx [ d A k ■ 



(2tt) 4 J J [k 2 + p 2 x(x - 1) + M 2 (l - x)] 2 

= _j^ nz-d/2) r 1 i 

(2ir) d T(2) J (1 - x) 2 - d / 2 (-p 2 x + M 2 ) 2 ~ d / 2 ' K> 

In the region of low momenta, where the chiral expansion is expected to converge quickly, the first factor in 
the denominator is of O(q ), while the second is of 0(q 2 ) so that the chiral dimension of the integral should 
be 0(q d ~ 2 ). We remark that at x = 1, the integrand develops a pole. Following we perform a change of 
variables 

II j2 

(10) 



A 4u 

p z 



so that 



J% 2 ) - T{2 - d/2) 1 du 1 (MY' 2 " 1 (ID 

Jq[P > (2n) d T(2) p 2 J dU (1 - au)2-rf/2(l _ u f-d/2 W«) • ^> 

We see that a complication arises due to the infrared singular piece ~ l/{k — p) 2 in the spin-projected parts 
of the propagator. Specifically, due to the term ~ au 2 , the denominator brings some additional a contribution 
when performing the integral. Jq{p 2 ) turns out to behave as M d ~ 4 instead of the expected M d ~ 2 . If one 
evaluates this integral straightforwardly, one finds that it contains a logarithmic contribution that indeed has 
the correct chiral behavior. The disturbing contribution is due to the factor (1 — au) which stems from the pole 
at x = 1 of the infrared singular part of the propagator. In fact, such contributions should be contained in the 
regular part, since they have nothing to do with the chiral expansion. To formally achieve this separation, one 
rewrites the integral as (d = 4 — 2e) 

l/a -y r^/ a 

du- N9 _ rf/9/1 7^—777 = / du(l-eln(l-au)-eln(l-u)) 



o (1 - au) 2 -*/»(l - u)- -,- jo 



l/a r l/a 



1 



- - (W " 1) + / dU (1 _ u) 2~d/2 ■ ( 12 ) 



(I 



From the second integral one then obtains the irregular part by letting l/a — > oo. Therefore, the separation of 
the integral Jq(p 2 ) into the irregular and the regular terms reads, Jq(p 2 ) — Io{p 2 ) + Rq(p 2 ), 

^ 2 r(2-d/2) i r°° i (M 2 )d/ 2-i n ^ 
/o(p) - J2^r r(2) V 2 Jo du (i-u) 2 - d / 2 (MJ ' (13) 



R°( P 2 ) ^ r(2 - d/2) { I 1 dx 

K oiP) - (2?r)d r(2) i / dx 



1 -1 



3C: 



(1 -x)) 2 - d ' 2 

The irregular part can now be worked out using the methodology developed in we find 



,/r .: P 2 .v ■ Mh 2 -'I ' lN) 



CV) = -^f2A + ^In^l=-%. (15) 
1oit z m z J p z 



P 2 



4 



k 



> — ■ > — > 

P p-k 

Figure 1: Feynman diagram for the nucleon self-energy with an spin-3/2 intermediate state. Solid, double and 
dashed lines denote nucleons, deltas and pions, in order. 



Here, the piece ~ A contains a pole term, 



,d-4 



A = i^i^4-5 (1Il47r+ ^ +1) 



(16) 



with "/e the Euler-Mascheroni constant and we have set the running scale of dimensional regularization, A, to be 
the nucleon mass [§]• Furthermore, the pion tadpole contribution A ff is given by the well-known loop integral 



d d k 



i J (2tt) 2 Ml 



(17) 



Of course, the scale that naturally appears in loop integrals with deltas is toa, not to, as used in Eas.l jl5ll6|l . 
However, the differences due to this choice of scale can be absorbed in polynomial (regular) terms, so that 
one can work with one scale only in the coupled pion-nucleon-delta system. However, to assess the theoretical 
uncertainty at a given order, one can set A = toa in the delta loop diagrams. Ea. (|15f) also shows that the 
non-propagating parts of the projected propagator will be absorbed in tadpole contributions, as it is expected 
since such contributions can be understood as a very heavy delta shrunk to a point-like vertex. 

Armed with this prescription, we can now set up a consistent expansion in small external momenta, pion masses 
and the 7VA mass splitting A — toa — to, in complete analogy to the small scale expansion formulated for heavy 
spin-3/2 fields 5 . We therefore refrain from discussing here in any detail the corresponding power counting 
and the construction of the chiral effective Lagrangian (most of this is available in the literature) . As discussed 
in detail in jjH], special care has to be taken to fulfill the decoupling theorem of QCD, which states that leading 
chiral singularities can not be modified by resonance fields. From here on, all small parameters (quark masses, 
external momenta, nucleon-delta mass splitting) are collectively denoted as e. 

5. We first calculate the leading one- loop contribution to the nucleon self-energy with an intermediate delta 
without any reference to a special regularization and with the full propagator, Eq.QJ, see Fig.^ This is done 
for the following reasons. First, such a calculation allows one to study different regularization procedures. 
Second, one can explicitly study the role of the non-propagating spin-1/2 components. Third, one can also 
easier compare with existing calculations in the literature. We have 



Sa(p 2 ) 



d 4 k 



4F2 J (2^)4 
The general structure of this integral is 



Qt{p-k) 



1 



Ml - k 2 



k v 



(18) 



w) = c ^ / (^( fc ) + B ^y+ C ( k W 



(19) 

i j y^/i ) - \ j 

with C = 9g\/(4F%), and the functions A(k), B(k) and C(k) can be straightforwardly evaluated. After 
integration, the self-energy can be written as 



£a(p 2 ) = A+rfB = A + Bm 



to) B 



(20) 
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where the functions A and B follow after integration of Ea. ll9|) . The nucleon mass shift due to the delta 
intermediate state is (setting Z = —1/2 for simplicity) 

5m = (A + B m)\ p 2 =m 2 

(d-2)C ( 1 2 . 1 , A A. , A A. 2l , 2 , 

= ti — tt — 5- -— to a toAa - -m A m Aa + -t(toa + to)toAa + — (toa + to) to Ji(m ) 
(a — 1) m A \ 2a 2 4 2 



Af 2 

+M 2 to 2 (to a + to) J (m 2 ) ^ A A - toM 2 Ji(to 2 ) 



2,1 ,1 2 

m + — to a to + 2 TO A 



-M>J 1 (m 2 ) + M 2 A T 



mA + to + ^ m 



(21) 



in terms of the loop functions 

Mp 2 ) 



A; 



and wc have defined 



d d k 


1 






(27T) 2 (Ml 


— k 2 )(m\ 


"(P- 


fc) 2 ) 


d d k 








(2tt) 2 (Af 2 


— k 2 )(m\ 


"(P- 


fc) 2 ) 


d d fc 1 








(2tt) 2 to a - 


-k 2 ' 






A = TOA ~ 


TO . 







(22) 
(23) 



In IR, we have Aa = 0, and Jo (to 2 ) (Ji(to 2 )) is of order 0(e) (0(e 2 )), so that to lowest (third) order 



5m = 



{d-2)C 
(d-l)m 2 A 



— A (to + toa) 2 to 2 Ji(to 2 ) + Af 2 to 2 (to + toa) Jq(to 2 ) 



3.9 



16tt 2 F 2 



A A 2 - >/ 2 In + 2(A 2 - A/ 2 ) 3 / 2 In + 



A 



A 2 



-1 - 



A /4 



2 V3 



-A - AT; 



(24) 



where the last equation is valid for A 2 > AT 2 (as it is the case in nature). This result agrees with the one in 
|3U) . Note also that for a more compact notation and easier comparison with the existing literature, we have 
retained some higher order terms in the formula expressing the mass shift in terms of the loop functions Jo and 
Ji- 

As noted above, we have not yet separated the spin-3/2 from the spin-1/2 components. Since the spin-1/2 
components do not propagate, their contribution should be entirely absorbed in polynomial pieces (contact 
terms). Utilizing the projection operators given in Eq.|j3J), the spin-3/2 contribution to the self-energy takes the 
form 

3/2^ _ (d-2)C f d A k 2 , 2 -(]/+m A -k) 

{2ttY K 1 1 P > ' {k-p) 2 {M 2 -k 2 ){m 2 A -{k-p) 2 )' 



JA * ' d-1 

which leads to the nucleon mass-shift (we give here only the result employing IR) 



(25) 



5m z ' 2 



{d-2)G 
m 2 A 



M 2 m 2 ( Jo(m 2 )(m A + to) — mJi(m 2 ) ) + -m 2 ( A(toa + to) — Af 2 ) (toa + to) Ji(to 2 ) 



+ ^ ( -A(m A + to) + Ml ) toJi(to z ) - M 2 m 2 ( J^(m z )(m A + to) - J"(m 2 )m 



+ ^-(M 2 + to 2 )(to a + to) J°(to 2 ) - ^(M 2 + to 2 ) 2 J^to 2 ) 



in terms of the new loop functions Jq(p 2 ) defined in Eq.@ and 



Pm J lip 



1 f d d k 



i J (2tt) 2 {p - k) 2 (M 2 - k 2 ) 



(26) 



(27) 



G 



We also have the loop function relation Ji(p 2 ) = (p 2 + M%) Jq (p 2 )/(2p 2 ) + A 1T /(2p 2 ) (in dimensional regular- 
ization). Noting now that the loop functions J® (to 2 ) and Jj 3 (to 2 ) are of order 0(e 2 ) in the chiral expansion, 
the lowest (third) order result for the mass shift derived from Ea. (|26|l agrees with the one given in Eq.J23J but 
differs of course in the higher order terms. It is also important to stress that the spin-3/2 contribution to the 
mass shift is independent of Z , as it should be. 

a/ V) 



Next, we consider the spin-1/2 contribution to the self-energy, E A 
form (we again give for simplicity the result for Z = —1/2) 

3si 1 



The corresponding mass shift has the 



8m 1 / 2 = 



2F%m% 256 



Ml 



M 2 

21n— f (2(m A 



<f-(2m - 



rrr 



M 4 

m A ) o- 



. , M 2 

+- (to + 4to a ) + 4— |(2to 



toa) 



M, 



m 

4\ 1 



(28) 



We have thus shown explicitly that these contributions can be completely absorbed into polynomial terms 
appearing in the chiral expansion of the nucleon mass beyond leading one-loop order |32l Ej , 



rriN = to — 4ciM n 



32ttF 2 



*i Ml In : 



Mir 



-k 2 M 4 + 0{Ml) 



(29) 



where c\, k\ and k 2 are (combinations of) dimension two and four low-energy constants and too is the nucleon 
mass in the chiral limit. Furthermore, all Z-dependence is of course included in these polynomial terms, too. 

From the mass shift, one can directly derive the delta loop contribution to the sigma term using the Feynman- 
Hellmann theorem, 

a , . o dSm 
<^(0)=M 2 — . (30) 

To leading order one has (using the full or the spin 3/2 part of the propagator): 



'ttN 



(0) 



M 2 (to + toa) 



C m 2 J (m 2 ) + 



1 1 to 2 + 2m 2 A 

6 47T 2 TO 2 



A (to + toa) 



CM 2 
~4^~ 



A In 



^ + v /A 2 -M^ln(A + . 



(31) 



Eq. H31[) naturally agrees with the result of Ref. but differs by a polynomial contribution in A from the 
one of jl(J| due to the use of a d-independent propagator in that paper. Again, for easier comparison we have 
retained some higher order terms in the upper expression in Eq. I|31l) . This concludes the formalism and we 
now turn to a numerical evaluation of the mass shift and the sigma term. 



6. We are now in the position to evaluate the delta loop contribution to the nucleon mass shift and the sigma 
term. We use g A = 1.267, M w = 139.57 MeV, F„ = 92.4 MeV and A = 271 McV. From Eq.© it follows 
that ca = 1.34. We remark that the resulting numbers scale linearly with the axial- vector coupling constant 
ca- The delta contribution to the nucleon mass shift and the sigma term are collected in Tabled We give 
the lowest (third) order result which is the same if one uses the full or the projected spin-3/2 propagator, as 
pointed out earlier. However, in IR the integrals also contain higher order pieces which can be retained (see 
e.g. the detailed discussion in H^D- In that case, one has to project onto the spin-3/2 pieces to get rid of 
the unphysical contribution from the spin-1/2 components. However, as shown in the table, the difference is 
numerically irrelevant for the observables considered here. Furthermore, these higher order corrections amount 
to a very small correction to the lowest order result. This, however, is not always the case, see e.g. ^S]. We 
also remark a certain sensitivity to the scale of dimensional regularization. This dependence, however, will be 
balanced by counter terms not considered here. The results shown in Tablen are n °t very different from earlier 
ones obtained in the heavy baryon approach |341 133j . which is expected since to this order neither the delta loop 
contribution to the nucleon mass shift nor to the sigma term is very sensitive to recoil corrections or similar 
effects. We stress again that a complete analysis of the nucleon mass or the sigma term would require a novel 
determination of certain low-energy constants to fulfill the requirements of decoupling. This, however, goes 
beyond the scope of this Letter. 



7 





5m [MeV] a£ N {0) [MeV] 


lowest order 

unexpanded, with full propagator 
unexpanded, with spin-3/2 propagator 


86.4 (108.1) -45.6 (-60.0) 
87.2 (112.7) -46.3 (-63.3) 
88.1 (113.7) -44.7 (-61.4) 



Table 1: Delta contribution to the nucleon mass shift and the sigma term. Lowest order refers to the third order 
IR result based on Eg. \2$ . Unexpanded means that the higher order terms from the various loop functions are 
retained. The propagator used is also given (in the lowest order case, using the full or the spin-3/2 part of the 
propagator leads to the same result, as discussed in the text). The numbers in the round brackets are obtained 
with A = rriA ■ 



7. In this Letter, we constructed a systematic infrared regularization for chiral effective field theories including 
spin-3/2 fields. To arrive at a consistent formulation, one has to project onto the spin-3/2 components of 
these fields. This projection leads to a new type of integrals as compared to the pure pion-nucleon theory. 
We have shown that the IR scheme of [5] can be extended systematically and furthermore, this procedure 
allows for a chiral invariant inclusion of external sources, like e.g. photons or weak currents. The contributions 
from the spin- 1/2 components get completely absorbed into the polynomial pieces of the effective Lagrangian. 
Furthermore, the frequently used off-shell parameters (X, Y, Z) are also contained in these polynoms and thus 
do not lead to observable consequences. It is now important to apply this method to observables, where Air 
loops are expected to play a significant role, like in all types of Compton scattering. Work along these lines is 
underway. 

TRH thanks the Universite Louis Pasteur for hospitality, where part of this work was done. 
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